Theory of high energy optical conductivity and the role of oxygens in manganites 
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Recent experimental study reveals the optical conductivity of Lai-^Ca^MnOa over a wide range 
of energy and the occurrence of spectral weight transfer as the system transforms from paramag- 
netic insulating to ferromagnetic metallic phase [Rusydi et ai, Phys. Rev. B 78, 125110 (2008)]. 
We propose a model and calculation within the Dynamical Mean Field Theory to explain this phe- 
nomenon. We find the role of oxygens in mediating the hopping of electrons between manganeses 
as the key that determines the structures of the optical conductivity. In addition, by parametrizing 
the hopping integrals through magnetization, our result suggests a possible scenario that explains 
the occurrence of spectral weight transfer, in which the ferromagnatic ordering increases the rate 
of electron transfer from 02 P orbitals to upper Mn es orbitals while simultaneously decreasing the 
rate of electron transfer from 02 P orbitals to lower Mn es orbitals, as temperature is varied across the 
ferromagnetic transition. With this scenario, our optical conductivity calculation shows very good 
quantitative agreement with the experimental data. 

PACS numbers: 



Introduction. Manganites have been the subject of ex- 
tensive studies since they have exhibited a wealth of fasci- 
taning phenomena such as the colossal magnetoresistancc 
(CMR), charge-, spin-, and orbital orderings, and transi- 
tion from paramagnetic insulator to ferromagnetic metal, 
as well as multifcrroic behavior [J-Q. Upon hole dop- 
ing, the transition from antiferromagnctic insulator to 
ferromagnetic metal has been argued to occur through a 
mixed-phase process [j| . Whereas for a fixed hole doping 
where ferromagnetic order is found, insulator to metal 
transition simultaneously occurs as temperature is low- 
ered across the ferromagnetic transition Q. R has been 
generally assumed and experimentally confirmed that the 
magnetic order in these systems is driven by the double- 
exchange interactions (o1.[7rll0|. However, explanation on 
other phenomena accompanying the ferromagnetic tran- 
sition seems to be far from complete, and remains as an 
open subject. 

Several theories on the insulator-metal (PM) transition 
accompa nying t he ferromagnetic transition have been 
proposed P^EDl- Although the details of the models and 
scenarios of the PM transition proposed by these theo- 
ries are quite different, they have similar idea suggesting 
that the Jahn- Teller distortion along with the electron- 
phonon interactions stabilize the insulating phase at high 
temperatures, which is broken by the ferromagnetic order 
below its transition temperatures. Phese theories, how- 
ever, have only addressed the static properties or low en- 
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ergy phenomena, as their models implicitly assume that 
low energy phenomena occuring in these materials are in- 
sensitive to possible high energy excitations. Many such 
models [lJ-[l9| typically consider only effective hoppings 
between Mn sites, while ignoring the electronic states in 
oxygen sites. On the other hand, models that included 
local interactions and hybridization in correlated mate- 
rials, might expect pronounced effects at higher energies 
that are connected to charge-transfer or Mott-Hubbard 
physics (2fJrf23j . Phus, the validity of such theories may 
have to be tested through experimental studies on the 
band structures and the optical properties over a wide 
range of energy. In that respect, experimental studies of 
optical conductivity of manganites as functions of tem- 
perature and doping in a much wider energy range be- 
come crucial. 

A recent study of optical conductivity by Rusydi et al. 
[24| . has revealed for the first time strong temperature 
and doping dependences in Lai-^Ca^MnOs for x = 0.3 
and 0.2. Phe occurrence of spectral weight transfer has 
been strikingly found between low (<3eV), medium (3- 
12eV), and high energies (>12eV) across PM transition. 
In fact, as the temperature is decreased, the spectral 
weight transfer appears more noticeably in the medium 
and high energy regions than it does in the low energy 
region. Observing how the spectral weight in each re- 
gion of energy simultaneously changes as temperature is 
decreased passing the ferromagnetic transition tempera- 
ture {Tfm)i one ma y suspect that there is an interplay 
between low, medium, and high energy charge transfers 
that may drive many phenomena occuring in manganites, 
including the PM transition. Phis conjecture is related 
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FIG. 1: (Color online) Simplified crystal structure of the 
model. The crystal structure is assummed cubic with the unit 
cell containing only one Mn and three O atoms. The O atoms 
are labeled 1,2, and 3 to distinguish the p— orbitals belonging 
to different O atoms used to construct our Hamiltonian. 



ther, with some additional argument, our calculation cap- 
tures qualitatively correctly the temperature dependence 
of the optical conductivity as the system transforms from 
paramagnetic to ferromagnetic phase. 

Model. As shown in Fig. [TJ we model the crystal struc- 
ture of Lai-zCa^MnOs such that each unit cell forms a 
cube with lattice constant a set equal to 1, and contains 
only one Mn and three O sites, thus ignoring the pres- 
ence of La and Ca atoms that we believe not to contribute 
much to the structures and temperature dependence of 
the optical conductivity. We choose 10 basis orbitals to 
construct our Hilbert space, which we order as the fol- 
lowing: |Mn e g U pper,t), |Mn e g i OW er,t), |Oi Pf), |0 2 pf), 
|0 3 pf), |Mn e g U pper,i), |Mn e g i owe .r,\), Oi p±), |0 2 p±), 
and |C<3 pi). Note that the distinction between | e g U p per ) 
and | e g i OW er) states is associated with the Jahn- Teller 
splitting. Using this set of bases we propose a Hamilto- 
nian: 



to the fact that the hopping of an electron from one Mn 
site to another Mn site can only occur through an O site. 

Considering the difference between the on-site energy 
of the manganese and that of the oxygen that could be 
about 5-8 eV [25j], the Mn-0 hoppings occur with high 
energy transfer. We hypothesize that if such high en- 
ergy hoppings can mediate a ferromagnetic order, then 
other low or high energy phenomena may possibly occur 
simulataneously. Thus, the mechanism of FM transition 
in the dc conductivity may not be completely separated 
from what appears as the decrease (increase) of the spec- 
tral weight in the medium (high) energy region of the 
optical conductivity, all of which together may be driven 
by the ferromagnetic ordering. Theories based on effec- 
tive low energy models which only consider Mn sites while 
ignoring O sites would not be able to address this. 

Motivated by the aforementioned conjecture, we de- 
velop a simple but more general model, in which oxygens 
are explicitly incorporated. In this paper, we propose 
our model and calculation of the optical conductivity of 
Lai_2;Ca x Mn03 within the Dynamical Mean Field The- 
ory, to explain the experimental results of Ref. [24[ . Our 
calculated optical conductivity shows that both oxygens 
and mangeneses play important roles in forming struc- 
tures similar to those of the experimental results. Fur- 



+ ^2 U u n Uit n ua + Y Uini if n ui 

ia icr 

~Y / JhS 1 .s 1 . (1) 



The first term in the Hamiltonian is the kinetic part, 
whereof r/^. is a row vector whose elements are the cre- 
ation operators associated with the 10 basis orbitals, and 
?7k is its hcrmitian conjugate containing the correspond- 
ing destruction operators. Here we consider that each 
Mn site contributes 4 e g orbitals (upper and lower each 
of which is with up and down spins), and the three O 
sites contribute 6 orbitals (3 from each site each of which 
is with up and down spins). [iJo(k)] is a 10x10 matrix 
in momentum space whose structure is arranged in four 
5x5 blocks corresponding to their spin directions as 
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where O is a zero matrix of size 5x5, and (referring to 
the choice of coordinates in Fig. [1]) 
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The diagonal elements of i?o(k)-|-m represent the local energies, while the off-diagonal elements represent the 
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hybridizations between orbitals. The first two diagonal 
elements, i.e. Ejt and —Ejt, correspond to the Mn 6j 
orbital energies which are split due to the presumedly 
static Jahn- Teller distortion. Each of the remaining three 
diagonal elements, i.e. E p , corresponds to the local en- 
ergy of the C>2p orbital. The parameter tj^-O (*M n -o) 
corresponds to hopping between the upper (lower) Mn ej 
orbital and the nearest C>2p orbital. Whereas to-o cor- 
responds to hopping between nearest C>2p orbitals. 

The second term in Eq. ([1]) represents the Coulomb 
repulsions between the upper and lower Mn e orbitals in 
a site. The third and forth terms represent the intra- 
orbital Coulomb repulsions. In this work, we take U u 
and Ui to be infinity, forbidding double occupancy in 
each of the lower and upper Mn e orbitals. Finally, the 
fifth term represents the double-exchange magnetic inter- 
actions between the local spins of Mn, S, formed by the 
strong Hund's coupling among three ti g electrons giving 
5=3/2, and the itinerant spins of the upper and lower 
Mn e electrons, s. Note that we use a well-accepted 
general assumption that the on-site Coulomb repulsion 
in each tig orbital and the Hund's coupling among the 
t2 g orbitals are so strong to keep the occupancy of the 
three tig levels fixed at high spin configuration. Thus 
the charge degrees of freedom of the three t2 g electrons 
become frozen, and the remaining degree of freedom to 
be considered is the orientation of the collective spin 3/2. 

Method. To solve our model, we use the Dynamical 
Mean Field Theory [26| . First, we define the Green func- 
tion of the system, which is a 10 x 10 matrix, 



[G(k, ; 



[ffo(k)] + [£(*)] 



(4) 



with z the frequency variable. Then, we coarse-grain it 
over the Brillouin zone as 



[G{i 
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In defining [G(k, z)], all the interaction parts of the 
Hamiltonian (all terms other than the kinetic part), are 
absorbed into a momentum-independent self energy ma- 
trix, [S(z)], which will be solved self-consistently. Note 
that in this algorithm, we need to go over the self- 
consistent loops in both Matsubara (z = iuj n + /i) and 
real frequency (z = uj + i0 + ). 

On taking U u and Ui to be infinity, to some ap- 
proximation, we forbid the double occupancies in states 
|Mn e g uppers), l Mn e g lower, \) , |Mn e g upper,]}, and 
|Mn e g lower. i) by throwing them out of our Hilbcrt space. 
To do this, according to the structure of Hamiltonian ma- 
trix in Eqs. ([2]) and ([3]), we multiply the weights of all 
the diagonal elements with indices 1, 2, 6, and 7, and all 
the corresponding off-diagonal elements connecting any 
pair of them by a half. Thus, after obtaining the ma- 
trix [G(z)] from Eq. ©, the effective [G(z)] (let's call 
it [G(z)] G ff) can be obtained by multiplying each of the 



following blocks of [G(z)] by a half, while keeping the 
remaining elements unchanged, that is 



Gn 


G12 


1 


Gn 


G12 


G21 


G22 


2 


G21 


G22 


Gi6 


G17 


1 


Gi6 


G17 


G26 


G27_ 


=> 2 


G26 


G27 


G6i 


G62 


1 


G6i 


G62 


G71 


G72 


=> 2 


G71 


G72 


G66 


G67 


1 


G&6 


G67 


G76 


G77 


=> 2 


G76 


G77_ 



(6) 



The "mean-field" Green function can then be extracted 
as 



[G(z)} = 



(7) 



Next, we construct the local self energy matrix, 
(z)], corresponding to the second and the fifth terms 
of the Hamiltonian. Here ni £ {0,1} is the occupation 
number of the lower Mn e orbital. The elements of the 
10 x 10 matrix [E ni (z)] are all zero except for the blocks 
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The local interacting Green function matrix is then cal- 
culated through 



[G ni (z)] = 



(9) 



where 6 and <f> are the corresponding angles representing 
the direction of S in the spherical coordinate. 

For each Mn site with a given m, the probability of 
Mn spin S having a direction with angle 8 with respect 
to the direction of magnetization (which is defined as the 
z— axis) is given by 
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d(cos0)e- s "< (cose) 
is the local partition function, and 

S ni (cos0) = -^lndct[G ni ( lW „)]e- 4U " 0+ (12) 



is the effective action. 



We need to average [G ni (z)] over all possible 9 and ni 
values as 



[G(z)] avo = (l-<n,» J d(cos9)P (cos9)[G (z)} 

+ (m) [ d(cos9)P 1 (cos9)[G 1 (z)}, (13) 



where (n{) is the average occupation of lower Mn e or- 
bital. The new self energy matrix is extracted through 



[£(*)] = [g(z)}- 1 - [G(z)j 



(14) 



Finally, we feed this new self energy matrix back into the 
definition of Green function in Eq. (U) , and the iteration 
process continues until [£(z)] converges. 

After the self-consistency is achieved, we can compute 
the density of states as 



DOS(u) = — ImTr[G(cj + i0+)] 



(15) 



We can also compute the optical conductivity tensor as 
^ f 'f(v,T)-f(v + u;,Ty 
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where [i'A(k)] = d[Ho(k)]/dk\ is the Cartesian compo- 
nent of the velocity matrix, [j4(k,i/)] = ([G(k, (jj + iQ + )} — 
[G(k,u — i0 + )])/(27ri) the spectral function matrix, and 
f(v,T) the Fermi distribution function. Note that the 
dimensional pre- factor ire 2 /(had) [26| with d = 3 is intro- 
duced to restore the proper physical unit, since the rest 
of the expression was derived by setting e = h = a = 1. 
In our model, the system is isotropic, and we are only 
interested in the tranverse components cr aa (oj) = <j(lj). 
which are equal for all a £ {x, y, z}. 

Results. Our calculated density of states is shown in 
Fig. [3J The parameter values used for this calculation are 



£jr=0.5eV, £ p =-6.5eV, i$,_ =1.2eV, 1 



(2) 
Mr, 



=0.8cV, 



t O -o=0.6eV, C/=10eV, Jjj=1.5eV, a=3.9345A, and T w 
194K (corresponding to /3 = l/T = 60 cV^ 1 ). These 
parameter values are chosen considering rough estimates 
given in other papers [HI, [l]| [25| and adjusted so as to 
give best agreement with the experimental optical con- 
ductivity data in Rcf. [24]. The DOS is normalized 
such that the integrated area is equal to 8, since in each 
unit cell there are 6 orbitals coming from oxygens and 
effectively 2 from manganeses, considering the restric- 
tion given by relation ([5]). The chemical potential is 
self-consistently adjusted to satisfy the electron filling of 
6+(l-a;)=6.7, mimicking the situation of Lai-^Ca^MnOa 
for x = 0.3. The structures of the density of states can 
be explained as the following. The three peaks labeled 
1,2,3 result from the fact that there are three oxygen 
atoms in a unit cell, where the degeneracy is broken into 
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FIG. 2: (Color online) Calculated density of states (DOS). 
See text for the parameter values used in this figure and the 
detailed explanation of the structures of the DOS. 



three levels by the hybridization between 2p orbitals of 
the neighboring oxygen atoms. The structures labeled 
with 4 through 9 result from the e g orbitals of mangene- 
ses. As shown in the figure, there are three mechanisms 
that split the Mn Cj states into 6 levels: static Jahn- Teller 
(JT) distortion, Coulomb repulsion U e ff) between lower 
and upper JT-split e g states, the double-exchange (DE) 
interaction between spins of electrons in the lower and 
upper JT-split states and the Mn spins formed by the 
Hund's coupling among the Mn t2 g electrons [28| . 

Figure [3] shows our calculated optical conductivity 
for T w 194K (> Tfm). The ferromagnetic transi- 
tion temperature for this set of parameters is roughly 
Tfm ~ 160K (based on extrapolation of the mean-field 
trend). The parameter values for T sa 194K are the same 
as those used in Fig. [2] In Figure [3j we demonstrate 
how we tune the the profile of the optical conductivity to 
achieve the best resemblance with the experimental data 
24 1 . It is important to note that our model is 



in Ref. 

not meant to address the dc conductivity, as we already 
anticipate that it cannot form an insulating (or nearly in- 
sulating) phase at T > Tfm , possibly due to not incorpo- 
rating electron- phonon interactions [HI, [ill [27j ■ Rather, 
our goal is to show how this simple model can capture 
qualitatively the general profile of the optical conductiv- 
ity from about 1 eV away from the Drude peak up to 22 
eV (the energy limit of the experimental data). 

On calculating the optical conductivity from Eq. (|16[) 
we introduce an imaginary self energy for the 02 P states, 
— iT, where r — l/T corresponds to the lifetime of the 
02 P states. The red curve in Fig. [3] shows the result if 
we use the self-consistent chemical potential, [i, in Eq. 
(fTB")) . Here, we observe that the resulting profile around 
the medium energy region (« 5-11 eV) does not satis- 
factorily resembles that of the experimental data in Ref. 
[24| . since some spectral weight seems to be missing in 
that region. We argue that the reason for this is related 
to the fact that our self-consistent chemical potential, 
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FIG. 3: (Color online) Calculated optical conductivity. Main 
panel: The red and blue curves represent the calculated opti- 
cal conductivites for different position of chemical potential, 
fi and fj.* , respectively, using only a small broadening (F=0.1 
eV). The black curve represents the result using chemical po- 
tential at n* with a bigger broadening (F=0.6eV). Inset: Re- 
gion in the density of states showing how the position of chem- 
ical potential is shifted. See text for detailed explanation. 



jj,, does not lie inside a pseudogap as it probably would 
if we incorporate electron-phonon interactions. In this 
model, wc only have a pseudogap that results from the 
double-exchange splitting, where /j, falls slightly to the 
right outside of this pseudogap. To remedy the missing 
of spectral weight, we slightly shift the position of chemi- 
cal potential to the left, i.e. from fi to /i* as shown in the 
inset of Fig. [3] Using this new chemical potential, /j,* , the 
resulting optical conductivity, shown in the blue curve, 
resembles the experimental data better. This suggests 
that the true chemical potential may actually lie inside 
a pseudogap similar to the situation as though it lies at 
fx*. (Note that, as long as considering the optical con- 
ductivity region about 1 eV away from the Drude peak, 
choosing fi* between 0.1 and 0.9 eV, i.e. around the val- 
ley, lead to similar results.) Although the profile of the 
blue curve is already better than the red one, it still has 
more pronounced stuctures than the actual experimental 
data does. To further tune the calculated optical conduc- 
tivity to better resemble the experimental data, we find 
that the overly pronounced structures can be broadened 
by enlarging the C>2 P imaginary self energy upto T = 0.6 
eV. The result after the broadening, which is shown in 
the black curve, looks very similar to the experimental 
results shown in Fig. 2(b) of Ref. (24| (replotted in the 
inset of Fig. 2]). This similarity in both magnitude and 
profile of the energy dependence may be a good measure 
of the validity of our model. 

Now we discuss how the model captures the spec- 
tral weight transfer when temperature is decreased from 
T > Tfm to T < Tfm- First, we divide the energy 
range into three regions: I (low:~l-3eV), II (medium:f«3- 
12eV), and III (high:>12eV), following the division made 



FIG. 4: (Color online) Spectral weight transfer in the optical 
conductivity. Main panel: Results of the model. Exclud- 
ing the region containing the Drude peak (0-leV), the energy 
range is divided into three regions I, II, and III. The black 
curve represents the optical conductivity in the paramagnetic 
phase, while the red, green, and blue curves correspond suc- 
cessively to lower temperatures in the ferromagnetic phase 
[29l |. The borders between regions I-II and II-III, denoted by 
the blue vertical dashed lines, are defined such that the curves 
are crossing at these energies. Inset: A replot of the corre- 
sponding experimental data from Ref. [24| for comparison. 



for the experimental data in Ref. [24[, except that we 
exclude the region around the Drude peak from our dis- 
cussion, since to obtain the correct values of conductivity 
in that region requires a more accurate description of the 
renormalized band structure around the chemical poten- 
tial. If we decrease temperature from the paramagnetic 
to ferromagnetic phase while keeping all the parameters 
constant, we find no significicant change in the optical 
conductivity, thus spectral weight transfer does not oc- 
cur in this way. If we inspect how Eq. (|16p determines 
the optical conductivity, we see that the change in op- 
tical conductivity may become more significant if either 
the spectral function, [j4(k, v)\, or the velocity operator, 
[u Q (k)], changes significantly while temperature changes. 
Within our model this can only be accomodated if we al- 
low some parameters to depend on temperature by some 
manner. By comparing the structures of optical con- 
ductivity and the corresponding DOS profile, it is clear 
that the spectral weight in the medium energy region 
comes mostly from transitions from C>2p to lower Mn e 
states, while in the high energy region from Oi v to up- 
per Mn e states. This fact may suggest that the hopping 



parameters tfy n _ and t^' depend on temperature. 
Furthermore, since the spectral weight transfer occurs 
most notably across and below Tfm, the temperature 
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dependence of t^n-O anc ^ ^Mn-O ma y ^ e reia ted to spin 
correlation. 

The actual interplay resulting in such a temperature 
dependence is believed to be very complicated, since 
it may involve orbital effects on the dynamic clcctron- 
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phonon coupling and spin correlation. In that regard, 
our present model, which is not an ab initio based 
model, cannot naturally capture these temperature ef- 
fects. Thus, to capture the plausible physics within 
our present model, we turn to the phenomenological ap- 
proach by parametrizing the totally non-trivial temper- 
ature effccs on hopping integrals through magnetization. 
In the simplest level, we may assume a linear dependence 
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of the hopping integrals t 
netization. Hence, we may write 
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4l-oW = 4l-o(0)(l + ci^) ! (17) 
<LM = 4l-o(0) (l + c 2 -J0 , (18) 



where M/M s is the ratio of magnetization to the satu- 
rated magnetization, and C\, and c 2 are constants. 
Using relations §Tj§ and (H]), taking i$ n _ o (0)=1.2 

eV, ^ n _ o (0)=0.8 eV, a 0.23, and c 2 « -0.35, at 
T=97K for which M/M s =0.357, for instance, we obtain 



that ^Mn-o ^ s enhanced to be 



1.3 eV, while t 
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suppressed to be « 0.7 eV. The results for four different 
temperatures are shown in Fig. 01 As shown in the main 
panel, our calculation shows that the spectral weight si- 
multaneously decreases (increases) in the medium (high) 
energy region of the optical conductivity as the system 
becomes ferromagnetic [30]. Our calculation also pro- 
duces a less noticeable decrease of the spectral weight in 
the low energy region as observed in the experimental 
data (see the inset). In both the main panel and the in- 
set, the black curve represents the optical conductivity in 
the paramagnetic phase, while the red, green, and blue 
curves correspond successively to lower temperatures in 
the ferromagnetic phase. If we define the positions of the 
borders between energy regions I-II and II-III such that 
all the curves are crossing at these energies, we obtain 
that theoretical values of these energies are similar to the 
experimental ones. Note that the temperatures varied in 
the theoretical and the experimental results should not 
be compared quantitatively, since the theoretical Tfm is 
about 100 K too small compared to the experimental one, 
possibly due to neglecting other possible exchange inter- 
actions in our model. Despite this, we believe that any 
improvement of Tpu by such additional terms would not 
change the physics presented in this paper. To show the 
difference in the density of states between paramagnetic 
and ferromagnetic phases, we display the spin dependent 
DOS for T w 194K and T « 97K in Fig. [5] 

To demonstrate further how the spectral-weight trans- 
fers in our model compare with the experimetal results, 
we display the relative spectral- weight changes for differ- 
ent regions of energy in Fig. [51 Comparing results in 
Fig. Eta) and (b), it is clear that for every region of en- 
ergy, I (low), II (medium), and III (high), (excluding 0-1 
eV), our calculations give exactly same trends as those 
shown by the experimental results. These suggest that 
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FIG. 5: (Color online) Spin-dependent density of states. The 
black and the red curves lie on top of each other as the spin- 
up and spin-down components of the DOS are identical in the 
paramagnetic phase. While, the green and the blue curves 
look quite distinct as DOS becomes polarized in the ferro- 
magnetic phase. 



the ingredients incorporated in our model are adequate 
to explain the occurrence of spectral-weight transfers in 
Lai-^Ca^MnOs in the energy range up to 22 eV. In that 
respect, one may argue, for instance, that the high-spin 
state (S = 3/2) of the ^g-electrons may become unstable 
as the system is optically excited by high-energy photons. 
Accordingly, transitions from high to low-spin states, or 
excitations of electrons from t2 g to e g levels may occur. 
Our present model does not incorporate those possibili- 
ties. However, our calculations prove that the model is 
capable to obtain the spectral- weight transfers with good 
qualitative agreement with the experimental results, thus 
suggesting that such other contributions may be minor 
or irrelevant. 

The inset of Fig. UJa) is to show that for 0-1 eV re- 
gion our result docs not agree with the experiment, since 
it does not capture the insulator-metal transition. As 
mentioned earlier, we argue that this is due to our model 
not incorporating the dynamic Jahn- Teller phonons and 
their interactions with electrons, which may be responsi- 
ble to form an insulating gap in the paramagnetic phase. 
The incorporation of such terms to improve our present 
model is under our on-going study. 

Conclusion. In conclusion, we have developed a model 
to explain the structures and the spectral weight transfer 
occuring in the optical conductivity of Lai-^Ca^MnOs 
for x = 0.3. The key that makes our model work in 
capturing the structures of the optical conductivity at 
medium and high energies is the inclusion of 2p orbitals 
into the model. 

Further, by parametrizing the hopping integrals 
through magnetization, our model captures the spectral 
weight transfer as temperature is decreased across the 
ferromagnetic transition temperature. Our calculation 
based on this phenomenological parameters suggests that 
the fcrromagnatic ordering increases the hopping param- 
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FIG. 6: (Color online) Relative spectral-weight changes. 
AW(T)/W from (a) our calculations, and (b) the experimen- 
tal results of Ref. [241] . for different regions of energy. AW/W 
is defined as spectral-weight difference AW = (cr{u), T) — 

a(u},TpM))dui normalized to W(Tpm) = J" 2 cr(aj, Tpm)(1uj, 
where in this case Tpm is 194 K in (a) and 310 K in (b). 
Positions of the theoretical and the experimental Tfm are in- 
dicated by vertical red dashed lines in each panel. Inset in 
(a) is comparison between AW (T) /W from the calculations 
(black filled diamonds) and from the experiments of Ref. [24T | 
(black empty diamonds) for 0-1 eV region. The horizontal 
green dashed line in the inset is just to highlight the zero 
position of AW/W. 



eter connecting the C>2p orbitals and the upper Mn eg or- 
bitals, while simultaneously decreasing the hopping pa- 
rameter connecting C>2p orbitals and the lower Mn es or- 
bitals. Although we have yet to check whether or not 
this scenario works in a more complete model incorporat- 
ing the dynamic electron-phonon coupling, we conjecture 
that this may be of important part that contributes to 
the mechanism of insulator to metal transition in man- 
ganitcs. 



Overall, our results demonstrate the strength of our 
model that one may have to consider as the minimum 
model before adding other ingredients in order to prop- 
erly explain the insulator-metal transition or other fea- 
tures in correlated electron systems such as manganites. 
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